Disordered Ground State and Magnetic Field-Induced Long-Range Order in an 
S=Z/2 Antiferromagnetic Honeycomb Lattice Compound Bi3Mn40i2(N03) 
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Bi3Mn40i2(N03), in which the Mn^^ ions carry 5=3/2, is the first honeycomb lattice system 
that shows no long-range magnetic order. Using neutron scattering, we have determined that short- 
range antiferromagnetic correlations develop at low temperatures. Applied magnetic fields induce a 
magnetic transition, in which the short-range order abruptly expands into a long-range order. 

PACS numbers: 75.25.-j, 75.30.Kz, 75.50.Ee 



Frustrated magnets show interesting phenomena origi- 
nating from the macroscopic ground state degeneracy. [l| 
Exotic states can be chosen as the ground state from 
many possible candidates. Usually, geometrically frus- 
trated systems consist of edge- and corner-shared trian- 
gles (triangular and kagome lattices, respectively) and 
tetrahedra (pyrochlore lattice). However, even unfrus- 
trated lattice systems, such as square or honeycomb lat- 
tice systems, can have frustrating interactions in the pres- 
ence of antiferromagnetic further-neighbor interactions. 

Magnetic field provides a perturbation to the degener- 
ate ground states of frustrated magnets, sometimes giv- 
ing rise to a novel phase transition. [2|-|j] The Cr-based 
spinels, consisting of the pyrochlore lattice, show a half- 
magnetization plateau that is driven by the strong spin- 
lattice coupling. 0, 0j 

The honeycomb lattice system with antiferromagnetic 
interactions has been previously studied. For example, 
SrL204 {L—Gd, Dy, Ho, Er Tm, and Yb) consists of dis- 
torted honeycomb latices. [6| SrEr204 shows an anoma- 
lous long-range magnetic order in which only one of 
the two crystallographically inequivalent Er sites orders, 
probably due to the combined effects of exchange inter- 
actions, dipolar interactions, and crystal- field. [7] 

Bi3Mn40i2(N03) has a trigonal structure (PS) with 
a=4.9692 A and c— 13.1627 A and consists of regular 
honeycomb lattices of the magnetic Mn^+ (nominally 
t2g;S = 3/2) ions, [sj Since a Jahn- Teller distortion is not 
observed in this compound, the orbital degrees of freedom 
appear to be unimportant. Therefore, this compound is a 
good model compound for an antiferromagnetic Heisen- 
berg honeycomb lattice system. The magnetic suscepti- 
bility shows a broad maximum centered around 70 K, a 
characteristic feature of a two-dimensional antiferromag- 
net. The Curie- Weiss temperature has a large value of 
©cw = —257 K, indicating strong antiferromagnetic in- 
teractions. Long-range magnetic order is not observed 
down to 0.4 K. This is much lower than |0cvf| and indi- 
cates the presence of strong frustration and suggests 



that the nearest-neighbor (Ji) and next-nearest-neighbor 
interactions (J2), as shown in Fig. [H compete in this 
compound. As far as we know, Bi3Mn40i2(N03) is the 
first honeycomb lattice system that shows no long-range 
order. 

We performed neutron scattering experiments in a 
powder sample of Bi3Mn40i2(N03). In zero magnetic 
field, it was confirmed that the compound does not show 
long-range magnetic order down to 3 K, indicating that 
the ground state is disordered. While no long range order 
is observed, short-range antiferromagnetic order, which 
is glassy in nature, develops at low temperatures. It 
was also found that the coupling between the two ad- 
jacent honeycomb planes that form a bilayer structure 
(Jc), as shown in Fig. [U is not negligible, showing 
that Bi3Mn40i2(N03) is a model compound for the bi- 
layer honeycomb lattice system. Under applied magnetic 
fields, this compound shows a commensurate long-range 
antiferromagnetic order, which corresponds to an expan- 
sion of the short-range order in zero field. This indicates 
that magnetic field perturbs the degenerate ground state 
and gives rise to a Neel order. 

A powder sample (4 g) of Bi3Mn40i2(N03) was pre- 
pared by hydrothermal synthesis. Q Although the sam- 
ple contains 6.7 weight % of Mn02 as shown in Ref. 
[a| , magnetic signals originating from the impurity phase 
were successfully removed. Elastic neutron scattering 
experiments were carried out on the thermal and cold 
neutron triple-axis spectrometers TAS-2 and LTAS, re- 
spectively, installed at the guide hall of JRR-3 at Japan 
Atomic Energy Agency. The horizontal collimator se- 
quences were guide-open-S-80'-open and guide-80'-S-80'- 
open on TAS-2 and LTAS, respectively. The fixed inci- 
dent neutron energies were 14.7 and 3.5 meV with the 
energy resolutions (AE) of 1.7 and 0.2 meV on TAS- 
2 and LTAS, respectively. Contamination from higher- 
order beams was effectively eliminated using PG and Be 
filters on TAS-2 and LTAS, respectively. A vertical mag- 
netic field up to 10 T was applied using a new type of 
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FIG. 1. (Color online) (a) Neutron powder diffraction pat- 
terns from Bi3Mn40i2(N03) in zero magnetic field at 4.8 
and 25 K. The broad peak originates from short-range an- 
tiferromagnetic correlations. The solid curves are the results 
of model calculations with (,ab=8 A and ^c~6 A. The thick 
horizontal bar represents the instrumental resolution, (b) The 
temperature dependence of the broad elastic magnetic signal 
measured with two different energy resolutions (Ai?=0.2 and 
1.7 meV). The two sets of the signals are normalized at the 
lowest temperature of each data. For all of the data shown 
here, the background data measured at ~50 K was subtracted, 
(c) Magnetic structure in the magnetic field-induced phase. 
Mn''"'" ions in 2x2x2 unit cells are shown. The magnetic mo- 
ments lie in the ab plane, although the direction of the mag- 
netic moments in the plane cannot be determined uniquely. 
Magnetic interactions Ji, J2, and Jc are also shown. 



split-pair superconducting magnet that uses cryocoolers. 
For the experiments in magnetic fields, the powder sam- 
ple was pressed into pellets to prevent the grains from 
orienting. Inelastic neutron scattering experiments were 
carried out on TAS-2. The horizontal collimator sequence 
was guide-80'-S-80'-40'. The fixed final neutron energy 
was 14.7 meV with an energy resolution of meV at 
E=0 meV. 

Figure [ija) shows neutron powder diffraction patterns 
for Bi3Mn40i2(N03) at 4.8 and 25 K. A peak much 
broader than instrumental resolution {AQ ~0.04 A~^) 
was observed at low temperatures. If the broad peak 
originates from a single honeycomb layer, an asymmet- 
ric peak is expected around (1, 0, 0) with a long tail at 
higher Q. This is not observed. Furthermore, there is no 
intense scattering at (1, 0, 2), indicating that the corre- 
lation between honeycomb bilayers is almost negligible. 
Since the observed peak position corresponds to (1, 0, 
1), only one bilayer unit with antiferromagnetic arrange- 
ment in a honeycomb layer and also between two honey- 
comb layers contributes to the scattering. The broad 
peak was simulated using the powder averaged three- 
dimensional (3D) squared Lorentzian with an isotropic 



FIG. 2. (Color online) Inelastic neutron scattering spectra in 
zero magnetic field in Bi3Mn40i2(N03) powder, (a) Image 
plot of the neutron scattering intensity in E-Q space at 5.5 
K. (b) Temperature dependence of the constant-energy scans 
measured at 2 meV. The solid curves are the results of the 
calculations with 3D squared Lorentzian. 



in-plane width, A/[{Q - Qwi)^ + 2£,J + where 
A, £^ab, and £,c are constant, effective correlation lengths 
in the ab plane and along the c axis, respectively. The 
powder averaging was performed based on the method as 
described in Ref. Q. Magnetic form factor was also in- 
cluded in the calculation. Since the peak width is much 
larger than the instrumental resolution, the resolution 
correction was not performed. The solid curves are the 
results of the calculations with ^q6=8 A and ^c=6 A. S^ab 
is larger than one hexagon and i^c is slightly longer than 
the distance between two honeycomb layers in one bi- 
layer. The latter result is consistent with the above result 
that the (1, 0, 2) peak is almost absent. It is noted that 
the spectra at 4.5 and 20 K can be reproduced equally 
well by the same length parameters, indicating that the 
correlation length saturates at low temperatures. 

Figure [TJb) shows the temperature dependence of the 
broad magnetic signal around the peak position. The 
elastic magnetic signal increases almost monotonically 
below the characteristic temperature (If), correspond- 
ing to a freezing temperature measured with the neutron 
energy scale. Tf strongly depends on the energy reso- 
lution of the neutrons. Tf becomes higher for measure- 
ments with broader energy resolution, with which signals 
from faster fluctuations are integrated. This is a typical 
behavior of the glassy magnetic systems in which spin 
fluctuations slow down gradually with lowering tempera- 
ture. This is consistent with the results of susceptibility 
measurements that show a hysteresis between measure- 
ments with field-cool (FC) and zero-field-cool (ZFC) pro- 
tocols. Ei 
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FIG. 3. (Color online) (a) Magnetization versus magnetic 
field in dc (T <10 K) and pulsed fields (T >15 K) for 
Bi3Mn40i2(N03) in a temperature range between 2 and 200 
K. (b) Temperature-magnetic field phase diagram. Filled and 
open circles show the data points determined from the mag- 
netization and neutron scattering measurements, respectively. 
LRO and SRO stand for the long-range and short-range or- 
dered phases, respectively. 



It was found that the ground state of Bi3Mn40i2(N03) 
is disordered. Our results indicate that Ji , J2, and Jc are 
dominant and there exist competing interactions between 
them. It is deduced from the short-range magnetic corre- 
lation that Ji and Jc should be antiferromagnetic. In ad- 
dition, J2 should be antiferromagnetic and compete with 
Ji . In the honeycomb lattice with Ji and J2 , the ground 
state is disordered when a(=J2/Ji)>l/6 and ~0.15 for 
classical spins and 5=3/2, respectively. [lll-[l3| 

Figure EJa) shows an image plot of inelastic neutron 
spectrum at r=5.5 K and H=Q T. Although there is 
no long-range magnetic order, relatively sharp magnetic 
excitations were observed. The spin-wave-like excitation 
rises from Q ~1.5 A~^, corresponding to the (1, 0, 1) 
position, and the band width of the excitation is ~9 
meV. An excitation gap originating from the magnetic 
anisotropy was not observed even with high energy res- 
olution experiments with AE ~0.2 meV, indicating that 
the interactions are isotropic. 

In order to clarify the temperature dependence of the 
spin fluctuations in more detail, temperature dependence 
of the constant-energy scans was measured at 2 meV as 
shown in Fig. (Hb). The peak width becomes broader 
with increasing temperature. The constant background 
intensity, which probably originates from paramagnetic 
scattering and a tail of nuclear incoherent scattering, in- 
creases with increasing temperature. The inverse peak 
width corresponds to the correlation length. The solid 
curve at 5.5 K is the result of calculations with S,ab=8 
A and ^c=6 A as in the elastic data shown in Fig. [Ija). 
The inelastic peak is broader than that of the elastic data, 
probably because the excitations are dispersive and also 
because the Q resolution is broader at 2 meV. Since the 
peak profile is symmetric at higher temperatures, indica- 
tive of isotropic spin correlations, the data at 50 and 100 
K were simulated using the powder averaged isotropic 3D 
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FIG. 4. (Color online) (a) Neutron powder diffraction pat- 
terns in Bi3Mn40i2(N03) with AE ~1.7 meV at H=0 and 
10 T and at T=3 K. At 10 T the broad magnetic peak is re- 
duced and sharp magnetic peaks appear. The magnetic field 
(b) and temperature dependences (c) of the elastic magnetic 
signal at Q=1.60 and 1.74 A"\ The signals at 1.60 and 1.74 
mostly originate from the short-range order and the 3D 
long-range order, respectively. In all the data shown here, the 
background signal was subtracted. 



squared Lorentzian. The solid curves are the results of 
the calculations with ^=3.3 and 2.5 A at 50 and 100 K, 
respectively. Even though the correlation length is very 
short, the peak is almost symmetric, suggesting that Jc 
is as large as Ji. 

From our inelastic neutron scattering experiments, the 
interactions can be estimated roughly. The maximum 
of the magnetic excitation, ~9 meV, corresponds to 
the sum of all the interactions multiplied by number of 
bonds and S'(=3/2), S'(3Ji -I- 6J2 + Jc), indicating that 
3Ji -I- 6J2 + Jc ^6 meV. Suppose that J2/Ji=0.15 corre- 
sponding to the critical value and Jc=0.5Ji, it would be 
estimated that Ji=2 Jc=1.4 meV and J2=0.20 meV. For 
comparison, Ji is estimated to be 3.0 meV from @cw 
{— —257 K), when J2 and Jc are not considered. The 
agreement is reasonably good, since Ji is reduced when 
J2 and Jc are included. 

Another characteristic feature in Bi3Mn40i2(N03) is 
that magnetization shows a small jump of ~0.1/iB/f.u. 
around 6-13 T below 20 K, as shown in Fig. [3)^ a), indi- 
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eating a transition to a metamagnetic phase. In order to 
elarify the nature of the magnetic field-induced phase, we 
performed neutron scattering experiments under applied 
magnetic fields. 

Figure SJa) shows the neutron powder diffraction pat- 
terns in Bi3Mn40i2(N03) at H=0 and 10 T. The spec- 
trum at T shows a broad peak originating from 
the short-range antiferromagnetic order as described 
above. With application of magnetic field at T—3 K, 
the broad magnetic intensity is reduced and resolution- 
limited sharp magnetic Bragg peaks appear above H ^6 
T, as shown in Fig. Sfb). This result indicates that long- 
range magnetic order is induced in magnetic field. Fig- 
urelH^c) shows the temperature dependence of the elastic 
magnetic signals at Q=1.60 and 1.74 at H^IO T. 
The signal originating from the long-range order vanishes 
and that from the short-range order recovers around 18 
K. 

Since the most intense peak is the (1, 0, 1) reflec- 
tion, the overall spin correlations do not change even in 
the magnetic field-induced phase. In the ab plane, the 
nearest-neighbor spins align antiparallel. The adjacent 
spins along the c axis also align antiparallel. Therefore, 
there is no frustration in the magnetic structure. Since 
the width of the magnetic field-induced sharp peak is al- 
most the same as the instrumental resolution, it is evident 
that the magnetic structure is long-ranged and commen- 
surate. From the results of the magnetic structural anal- 
ysis, we determined that the spins lie in the ab plane and 
that the ordered moment is '^1.8/iB at 3 K. The magnetic 
structure is shown in Fig. [Ijc). The moment size is about 
2/3 of the full moment 3/iB for the Mn^"*" ion as reported 
in Ref. Q. As shown in Fig. Ufa), there still exists a 
broad signal around 1.6 at 10 T. This is the remanent 
of the broad peak observed in the zero magnetic field. 
This indicates that the field-induced transition depends 
on magnetic field direction and that only the magnetic 
moments perpendicular to the magnetic field contribute 
to the ordering and the others remain disordered. It is 
puzzling why the transition depends on the field direc- 
tion, although the compound is an isotropic Heisenberg 
system. The Dyaloshinski-Moriya (DM) interaction may 
be an origin of the anisotropy, since the crystal structure 
(P3) does not have inversion symmetry. 

The magnetization measurements show that the meta- 
magnetic phase has a ferromagnetic component of 
^O.l/ie/f.u. A canting of spins probably occurs in the 
field-induced phase, although the small canted compo- 
nent cannot be detected in the present neutron diffraction 
measurements. It is possible that the canting is driven by 
the DM interaction, as described above. The stabiliza- 
tion of the magnetic ordering due to the spin canting in 
the magnetic moment is probably the origin of the field 



induced long-range ordering. 

The magnetic field-temperature phase diagram are 
summarized in Fig. ^h). We observed the magnetic 
field-induced phase transition, in which the short-range 
order expands into the long-range Neel order by apply- 
ing magnetic field. This indicates that this compound 
is located near the boundary of the ordered and disor- 
dered phases with a close to 0.15. It was reported in 
pyroclore systems Tb2Ti207 and Yb2Ti207 that mag- 
netic field along a particular direction induces a mag- 
netic transition to a long-range ordered phase although 
the detailed structure is unknown. [^, The transition 
temperature increases with increasing magnetic field in 
these compounds. As shown in Fig. [3Kb), the same be- 
havior was observed in Bi3Mn40i2(N03). This indicates 
that the collinear Neel state becomes more stabilized at 
higher temperatures. This may suggest that thermal fluc- 
tuations stabilize the Neel order, where "order by disor- 
der" mechanism is at work. 

We have observed disordered ground state and 
the magnetic field-induced transition to commensu- 
rate Neel order in the bilayer honeycomb system 
Bi3Mn40i2(N03). We emphasize that the zero-field and 
field-induced spin correlations in the honeycomb lattice 
system were studied experimentally for the first time. 
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